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Abstract. We revisit several maximization problems for geometric networks design under 
the non-crossing constraint, first studied by Alon, Rajagopalan and Suri (ACM Symposium 
on Computational Geometry, 1993) . Given a set of n points in the plane in general position 
(no three points coUinear), compute a longest non-crossing configuration composed of 
straight line segments that is: (a) a matching (b) a Hamiltonian path (c) a spanning tree. 
Here we obtain new results for (b) and (c), as well as for the Hamiltonian cycle problem: 

(i) For the longest non-crossing Hamiltonian path problem, we give an approximation 
algorithm with ratio « 0.4829. The previous best ratio, due to Alon et al., was 
I/tt ~ 0.3183. Moreover, the ratio of our algorithm is close to 2/n on a relatively broad 
class of instances: for point sets whose perimeter (or diameter) is much shorter than the 
maximum length matching. The algorithm runs in OirJ^'^ logn) time. 

(ii) For the longest non-crossing spanning tree problem, we give an approximation 
algorithm with ratio 0.502 which runs in O(nlogn) time. The previous ratio, 1/2, due to 
Alon et al., was achieved by a quadratic time algorithm. Along the way, we first re-derive 
the result of Alon et al. with a faster 0(n log n)-time algorithm and a very simple analysis. 

(iii) For the longest non-crossing Hamiltonian cycle problem, we give an approximation 
algorithm whose ratio is close to 2/7r on a relatively broad class of instances: for point sets 
with the product ( diameter x convex hull size ) much smaller than the maximum length 
matching. The algorithm runs in 0(n^'''' log n) time. No previous approximation results 
were known for this problem. 



1. Introduction 

Self-crossing in planar configurations is typically an undesirable attribute. Many struc- 
tures studied in computational geometry, in particular those involving a minimization con- 
dition, have the non-crossing attribute for free, for instance minimum spanning trees, min- 
imum length matchings, Voronoi diagrams, etc. The non-crossing property usually follows 
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from the triangle inequality. Alon et al. [3] have considered the problems of computing (i) 
the longest non-crossing matching, (ii) the longest non-crossing Hamiltonian path and (iii) 
the longest non-crossing spanning tree, given n points in the plane. Although they were 
unable to prove it, they suspected that all these problems are A^P-hard. The survey articles 
by Eppstein [HI pp. 439] and Mitchell [IT, pp. 680] list these as open problems in the area 
of geometric network optimization. The problem of approximating the longest non-crossing 
Hamiltonian cycle is also of interest and wide open [U pp. 338]. 

Without the non-crossing condition explicitly enforced, the problem of minimizing or 
maximizing the length of a spanning tree, Hamiltonian cycle or path, perfect matching, tri- 
angulation, etc. has a rich history. However if such structures are required to be non-crossing 
much less is known, in particular for the maximization variants. While for minimization 
problems, the non-crossing property comes usually for free via the triangle inequality, in 
contrast, for maximization problems, the non-crossing property conflicts directly with the 
length maximizing objective. This is another reason why these problems are interesting to 
study. 

Related work. The existence of non-crossing Hamiltonian paths and cycles in geometric 
graphs has been studied in pi [5]. Various Ramsey-type results for non-crossing spanning 
trees, paths and cycles have been obtained in [11] and jl2l . The Euclidean MAX TSP, the 
problem of computing a longest straight-line tour of a set of points, has been proven NP- 
hard in dimensions three or higher [9], while its complexity in the Euclidean plane remains 
open [14j. In contrast, the shortest non-crossing matching and the shortest non-crossing 
spanning tree are both computable in polynomial time jU [14], as they coincide with the 
shortest matching and the shortest spanning tree respectively. 

Definitions and notations. A set S of points in the plane is said to be in general 
position if no three points are collinear. General position will be assumed throughout this 
paper. Given a set of n points in the plane, the results of Alon al. are as follows: (i) A 
non-crossing matching whose total length is at least 2/7r of the longest (possibly crossing) 
matching can be computed in 0(n'^/^ log n) time, (ii) A non-crossing Hamiltonian path 
whose total length is at least l/vr of the longest (possibly crossing) Hamiltonian path can 
be computed in 0{rJ/'^ logn) time, (iii) A non-crossing spanning tree whose total length is 
at least n/{2n — 2) > 1/2 of the longest (possibly crossing) spanning tree can be computed 
in 0{'n?) time. The running times have been adjusted to reflect the current best upper 
bound of 0(n^/^) on the number of halving lines as established by Dey [6]. 

A geometric graph G is a pair {V, E) where y is a finite set of points in general position 
in the plane, and E' is a set set of segments (edges) connecting points in V . The length of 
G, denoted L{G), is the sum of the Euclidean lengths of all edges in G. The graph G is said 
to be non-crossing if its edges have pairwise disjoint interiors (collinear triples of points are 
forbidden in order to avoid overlapping collinear edges). 

For a point set S, let conv(S') be the convex hull of S, and let P = P{S) denote the 
perimeter of conv(S'). Denote hy D = D{S) the diameter of S and write n = \S\. Let 
MopT be a longest (possibly crossing) matching of S, and let MQprp be a longest non- 
crossing matching of S; observe that for odd n, Mqpt is a nearly perfect matching, with 
(n — l)/2 edges. Let -ffoPT be a longest (possibly crossing) Hamiltonian path of S, and let 
-^OPT ^ longest non-crossing Hamiltonian path of S. Let Toft be a longest (possibly 
crossing) spanning tree of 5, and let TQprp be a longest non-crossing spanning tree of S. 
Finally, let QoPT be a longest (possibly crossing) Hamiltonian cycle of 5, and let Qqpt 
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be a longest non-crossing Hamiltonian cycle of S. The following inequalities are obvious: 
L(Mopt) < L{Hopt) < L(roPT). 

Given a set 5 of n points in the plane, a line i going through two points of S is called a 
halving line if there are [(n — 2)/2j points on one side and \{n — 2)/2] points on the other 
side |T3]. A bisecting line £ of 5" is any line that partitions the point set evenly, i. e., neither 
of the two open halfplanes defined by I contains more than n/2 points of 5 [7]. Observe that 
any halving line of S is also a bisecting line of 5. Any bisecting line of 5 yields (perhaps 
non-uniquely) a bipartition S = RU B, with i? n -B = 0, — \B\\ < 1, with R contained 
in one of the closed halfplanes determined by i, and B contained in the other. We call 
S = RUB a. linearly separable bipartition, or balanced partition of S. Observe that for any 
non-zero direction vector v, there is a bisecting line orthogonal to v, see [71 Lemma 4.4]. 
Two bisecting lines are called equivalent if they can yield the same balanced partition of 
S. It is well known that the number of non-equivalent bisecting lines of a set is of the 
same order as the number of halving lines of the set, and any balanced bipartition can be 
obtained from a halving line [3, pp. 67]. 

Our results are summarized in the following three theorem^. 

Theorem 1.1. (i) For the longest non-crossing Hamiltonian path problem, there is an 
approximation algorithm with ratio ~ 0.4829 that runs in 0{n'^^^ logn) time. 

(ii) Given a set of n points in the plane, one can compute a non-crossing Hamiltonian path 
H in 0{n'^/^\ogn) time such that L{H) > ^L{Hopt) ~ ^- -^^ particular, if the point set 
satisfies the condition — < 6L{Hopt) for some small 6>0, then L{H) > (- - 5)L{Hopt)- 

(iii) Alternatively, one can compute a non-crossing Hamiltonian path H in 0{nlogn/ y/e) 
time, such that L{H) > (1 - e)^L{HoPT) - ^• 

Theorem 1.2. For the longest non-crossing spanning tree problem for a given set of n 
points in the plane, there is an approximation algorithm with ratio 0.502 and 0(n log n) 
running time. More precisely, the algorithm computes a non-crossing spanning tree T such 
that L{T) > 0.502 • L{Topt)- 

Theorem 1.3. Given a set S of n points in the plane, with |conv(S')| = h: 

(i) One can compute a non-crossing Hamiltonian cycle Q in 0(n'^/'^ log n) time such that 
L{Q) > ^L(Qopt) — (2/i — 1)^- In- particular, if the point set satisfies the condition 
{2h - l)f < 5L{Qopt) for some small 5 > 0, then L{Q) > (f - 5) L{Qopt)- 

(ii) Alternatively, one can compute a non-crossing Hamiltonian cycle Q in 0(n^ log n) time 
such that L{Q) > ^L{Qopt) - {h + 2)|. 

(iii) Alternatively, one can compute a non-crossing Hamiltonian cycle Q in 0(n log n/y^) 
time, such that L{Q) > (1 - £)^L{Qopt) - (2/i - l)f • 

2. The Hamiltonian path 

In this section we prove Theorem 11.11 Let S = {pi, . . . ,Pn}- We follow an approach 
similar to that of Alon et al. using projections and an averaging argument, in conjunction 
with a result on bipartite embeddings of spanning paths in the plane. Abellanas et al. [U 
Theorem 3.1] showed that every linearly separable bipartition S = RUB with | — | < 1, 
admits an alternating non-crossing spanning path such that the edges cross any separating 
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line i at points ordered monotonically along i. Such a Hamiltonian path can be computed 
in 0(n log n) time. Their algorithm computes the same Hamiltonian path for any two 
equivalent halving lines, that is, the alternating path depends on the bipartition only rather 
than the separating line. 

We now recall the algorithm of Abellanas et al. [1]; see Fig. H] for an example. Let 
S = RU B with ||i?| — \B\\ < 1 be the red-blue bipartition given by a vertical line £: R on 
the left, B on the right. Their algorithm constructs an alternating path A in the following 
way: Let rb be the top red-blue edge of the convex hull conv(5), called the top bridge. If 
|it;| > \B\, set A := {r}, if \R\ < \B\, set A := {b}, else set A to {r} or {b} arbitrarily. 
At every step, recompute the top bridge rb of S \ A, and add r to j4 if the last point in A 
was blue, or add 6 to A if the last point in A was red. As pointed out by the authors, the 
resulting path A is non-crossing because A is disjoint from the convex hull of 5" \ yl at each 
step. 

We improve the lower bound of Alon et al. by computing the longest Hamiltonian path 
corresponding to a bipartition and a Hamiltonian path of length at least the perimeter of 
the convex hull, and returning the longest of the two. 

Lemma 2.1. For a point set S, \S\ = n > 31, a non-crossing Hamiltonian path of length 
at least P{S) can be computed in O(nlogn) time. The bound on the length is best possible. 

Consider a geometric graph G = {V, E), and a point g ^ F, so that VU{q} is in general 
position. We say that q sees a vertex v \i the segment qv does not intersect any edge 
of G. Similarly, we say that q sees an edge e £ E, if the triangle formed by v and e does 
not intersect any other edge of G. We make use of the fact that if n is even then the two 
endpoints of an alternating path are on opposite sides of the separating line £. If n is odd, 
we first construct an alternating path for a specific subset of n — 1 points, and then augment 
it to a Hamiltonian path on all n points using the following lemma. 

Lemma 2.2. Let S = RU B with \ \R\ — \B\\ < 1, be a linearly separable bipartition given 
by line i. Let q G S, and A' be a non-crossing alternating path on S \ {q} such that its 
(consecutive) edges cross i at points ordered monotonically along £. Then q sees one edge of 
A' and consequently, A' can be extended to a Hamiltonian path A on S, with L(A') < L(A). 
The path A can be computed in 0{n) time, given A' . 

Fix a Cartesian coordinate system P. Let k be the number of halving lines of S, denote 
the angles they make with the x-axis of T by < oi < . . . < vr. By relabeling the points 
assume that the optimal path is -ffoPT = Pi,P2, ■ ■ ■ ,Pn- For two points Pi,Pj £ S, let /3ij 
be the angle in [0, vr) formed by the line through pipj and the x-axis. If n is odd, then a 
bisecting line of direction a (for any a) must be incident to at least one point of S, and 
denote an arbitrary such point by g^. 
Algorithm Al: 

Step 1 . Compute a non-crossing Hamiltonian path Hi of length at least P{S) , by Lemma[2T] 
Step 2. If n is even, then for all non-equivalent bisections of S (i.e., for all balanced bi- 
partitions of S) , compute a non-crossing alternating path using the algorithm of Abellanas 
et al. [1], and let the longest such path be H2. If n is odd, then for all non-equivalent bi- 
sections of S, compute a non-crossing alternating path of the even point set S \ {qa} using 
the algorithm of [IJ and let the longest such path be H'2. Augment H2 with vertex q^ by 
Lemma |2. 21 to a Hamiltonian path H2. 
Step 3. Output the longest of the two paths Hi and H2. 



LONG NON-CROSSING CONFIGURATIONS IN THE PLANE 



315 



By Lemma l2.ll the running time of Step 1 is O(nlogn). Since the number of halving 
hnes of an n-element point set is 0(n^/^) and all can be generated within this time [6], the 
running time of Step 2 is 0{'n7^^ logn), consequently the total running time of Al is also 
OirJl'^ logn). 

We proceed with the analysis of the approximation ratio. For simplicity, we assume 
that n is even. The case of n odd is slightly different. For each a G [0,7r), let Fq, be a 
(rotated) coordinate system, obtained from F via a counterclockwise rotation by a, and 
with the y-axis dividing evenly the point set S. Let Xj be the x-coordinate of point pi with 
respect to Fq. For a given a, let be a non-crossing alternating path with respect to a 
balanced bipartition induced by the y-axis of Fq,, as computed by the algorithm. There are 
0(1) balanced bipartitions given by any halving line of S. Recall that does not depend 
continuously on a; it depends only on the discrete bipartition. However, the coordinates of 
the points depend continuously on a. Assume that Ha = Pa{i)jPa(2)i ■ ■ ■ iPu(n)-, where a is 
a permutation of [n]; here a depends on the bipartition (hence also on a). Let Wa denote 
the width of S in direction a, that is, the width of the smallest parallel strip of direction a 
that contains S. By projecting on the x-axis of Fq, we get 

n 

L{Ha) > \ + 21x^(2)1 + . . . + 2|x^(„_i)| + |x<^(„)| = 2^ - - \x^^n)\ 

i=l 

n—l n—1 
= + + kll + \Xn\ - - \Xa{n)\ > + \xj+l\) - Wa 

j=l 3=1 
n-1 

^ |PiPj+i|| cos(/?jj+i - a)\ - Wq (2.1) 
i=i 

In the 2nd line of the above chain of inequalities, we use the fact that Pa{i) and Pa(n) 
lie on opposite sides of i, since n is even, hence + < \Pa{i)Pa(n)\ ^ 

In the 3rd line, we make use of the following inequality: for any two points Pi,Pj E S, 
\piPj\ \ cos(/3jj — a)| < |xj| -|- with equality if and only if the two points lie on opposite 
sides of the y-axis of Fq. Recall: for even n, H2 is the longest of the 0{k) Hamiltonian 
non-crossing paths i^Q. over all 0{k) balanced bipartitions of S. (A given angle Ofj yields 
0(1) balanced partitions, and corresponding alternating paths denoted here Hai-) We thus 
have for each a G [0,7r): 

n-1 

L{H2) > ^ \pjPj+i\\cos{Pjj+i - a)\ - Wa- 
i=i 

Note that 



/"TT /"TT 

/ I cos(/3jj_(.i — a)\ da = \ cos a\ da = 2, 
Jo Jo 



and according to Cauchy's surface area formula, we have Jq Wa da = P{S). By integrating 
both sides of the previous inequality over the a- interval [0, vr], we obtain 

n— 1 

7rL(i/2) > 2 \PjPj+i\ - PiS) = 2L{Hopt) - P{S), 
i=i 

HH2) > -L{Hovt) - (2.2) 

TT TT 
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We now improve the old approximation ratio of ^ ~ 0.3183 to ^jq^j ~ 0.4829, by balancing 
the lengths of the two paths computed in Step 1 and Step 2. Set c = 

Case 1: L(ffoPT) < cP{S). By considering the path computed in Step 1, we get a 
ratio of at least 

L{H,) ^ PjS) ^ P{S) _ 2 

L{Hopt) - L{Hopt) - cP{S) TT + 1 • 

Case 2: -L(i7oPT) ^ cP{S). By considering the path computed in Step 2 (inequality 
(j2.2p ). we get a ratio of at least 

HH2) ^ |L(goPT) - ^P{S) >2_^^2/_ ^\ ^ _J_ 
L{Hopt) ~ L^Hopt) ~7r ctt tt \ vr + ly vr + l 

Observe that if the point set satisfies the condition < SL^Hqpt:), then by (|2.2p . 
we have 

L{H) > -L{Hopt) - dL{HoPT) =(--s] L{Hopt). 

TT J 

This concludes the proofs of parts (i) and (ii) of Theorem II. 1[ 

(iii) With the same approach as in [3j, a Hamiltonian path of length at least (1 — 
£)^L{Hopt) — can be found by considering only b/^/e angles 6i = for i = 

0,1,..., [b/^/e\, where 6 is a suitable absolute constant. The resulting running time is 
0{nlogn/^/e). This concludes the proof of Theorem ll.il 

3. The spanning tree 

In this section we prove Theorem II. 2 [ Let S = {pi, . . . where pi = {xi,yi). Given 
a point p G S, the star centered atp, denoted Sp, is the spanning tree on S whose edges join 
p to all the other points. Since S is in general position, Sp is non-crossing for any p £ S. 
An extended star centered at p is a spanning tree of S consisting of paths of length 1 or 2 
(edges) connecting p to all the other points. See Fig. [TJ While the star centered at a point 
is unique, there may be many extended stars centered at the same point, and some of them 
may be self-crossing. In particular Sp is also an extended star. 




Figure 1: A star (left) and a non-crossing extended star (right) on a same point set, both centered 
at the same point p. 

The algorithm of Alon et al. computes the n stars centered at each of the points, and 
then outputs the longest one. The algorithm takes quadratic time, and the analysis shows a 
ratio of 2n-2 (which tends to 1/2 in the limit). Their algorithm works in any metric space. 
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As pointed out by Alon et al., the ratio 1/2 is best possible (in the hmit) for this specific 
algorithm. We first re-establish the 1/2 approximation ratio using a faster algorithm, and 
also with a simpler analysis. Our algorithm works also in any metric space; however in this 
general setting, the running time remains quadratic. 

Algorithm A2: Compute a diameter of the point set, and output the longest of the two 
stars centered at one of its endpoints. 

Obviously the algorithm runs in O(nlogn) time, with bottleneck being the diameter 
computation [15]. Let ab be a diameter pair, and assume w.l.o.g. that \ab\ = 1. The ratio 
1/2 (or even 2-0-2 ) follows from the next lemma in conjunction with the obvious upper 
bound L(Topt) < n (or L(Topt) 1^ n — 1). 

Lemma 3.1. Let Sa and Sb be the stars centered at the points a and b, respectively. Then 
L{Sa) + L{St) > n. 

Proof. Assume that a = pi, b = p2- For each i = 3, . . . ,n, the triangle inequality for the 
triple a, 6, pi gives 

\api\ + \bpi\ > \ab\ = 1. 

By summing up we have 

n 

L{Sa) + L{Sb) = X^da^'^l + \bPi\) + 2|a6| > (n - 2) + 2 = n. _ 

j=3 

We now continue with the new algorithm that achieves a (provable) ^ + approxi- 
mation ratio within the same running time 0(n log n). 

Algorithm A3: Compute a diameter ab of the point set, and output the longest of the 5 
non-crossing structures Sa, Sb, Sh, Ea, Eb, described below. 

Assume w.l.o.g. that the ab is a horizontal unit segment, where a = (0, 0) and b = (1,0). 
Let h = {xh,yh) be a point in S with a largest value of \y\. By symmetry, we can assume 
that yh > 0. Sa, Sb, and Sh are the 3 stars centered at a, b, and h respectively. Ea, resp. 
Eb, are two non-crossing extended stars centered at a, resp, b; details to follow. Each of 
the five structures can be computed in 0(n log n) time, so the total execution time is also 
0{n log n). 

Set 6 = 0.05, w = 0.6, t = 0.6 and z = 0.48, and refer to Fig. [3 Let h, £2, 4, and £4, 
be four parallel vertical lines: ii : x = 0, £2 : x = 0.2, £3 : x = 0.8, £4 : x = 1. Obviously, 
all points in S lie in the strip bounded by £1 and £4. Let Vm be the vertical parallel strip 
symmetric about the midpoint of ab and of width w. We refer to Vm as the middle strip; Vm 
is bounded by the vertical lines £2 and £3. Let Va and Vb be the two vertical strips of width 
0.2 bounded by £1 and £2, and by £3 and £4 respectively. Let c = {x^, Vc) be the intersection 
point between £3 and the circular arc 7a of unit radius centered at a and sub-tending an 
angle of 60°. We have Xc = 0.8 and 

Vc = Vl - 0.82 = 0.6 = t. 

We now describe the two extended star structures Ea and Eb. See also Fig. [3] for an 
example. To construct Ea, first compute the order of visibility of the points in Vb from point 
a by sorting. Then connect a with each point in the right strip Vb. Note that b G Vb, thus 
Vb / 0. Call S'a the resulting star. The edges of this star together with the vertical line £3 
divide Va U Vm into convex regions (wedges with a common apex a) ordered top-down. The 
subset of points in each wedge can be computed using binary search in overall O(nlogn) 
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y = 0.6 



V,„ 



Figure 2: A diameter pair a, 6 at unit distance, and tire tliree vertical strips Va, Vm, and Vf,. The 
two circular arcs 7^ and 7f, of unit radius centered at a and h intersect at the point 
(1/2, \/3/2). All points of S above ab lie in the region bounded by a6, 7^ and 76. 



time (over all wedges). 5^ is extended (augmented) as follows. In each wedge, say pag, all 
points are connected either to a or to depending on the best (longest) overall connection 
cost. We denote the resulting extended star structure by Ea- The construction of is 
analogous. It is clear by construction that both Ea and E^ are non-crossing. 

Lemma 3.2. For each p £ S, let dmaxijp) denote the maximum distance from p to other 
points in S. Then 



L{Topt) < 



.i=l 



Proof. Consider Tqpt rooted at a and drawn as an abstract tree with the root at the top in 
the usual manner. Let tt(v) denote the parent of a (non-root) vertex v. Uniquely assign each 
edge 'ir{v)v of Topt to vertex v. Obviously, L{-k{v)v) < dmnxiv) holds for each edge in the 
tree. By adding up the above inequalities, and taking into account that (imax(a) = |a&| = Ij 
the lemma follows. ■ 

Lemma 3.3. Assume that Y17=i IVil ^ for some positive constant 5 < 1. Then 



L{Sa) + L{Sb)>2n\l^ + 5^. 



Lemma 3.4. Let Ua and denote the number of points in the left and right vertical strips 
Va and Vh- Then L{Ea) > ^-^{n + ni,), and similarly L{Ei,) > ^-^{n + Ua)- Consequently 
L(Ea) + L{Ei,) > ^-|^(2n + Ua + rif,). Ea and E^ can be constructed in O(nlogn) time. 

Proof. The distance between li and ^3 is ^-5^. By an argument similar to that in the proof 
of Lemma |3. 11 the connection cost for a wedge with m points is at least ^-^m. Therefore 
the total length of Ea is 

T / \ ^ + W 1 + W 1+W 

L{Ea) > — - — Ub H — [n - Ub) = — —{n + rib)- 

The estimation of L[Eb) is analogous. The running time has been established previously. ■ 
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Lemma 3.5. Assume that Y17=i 1^*1 — ^^'^ Hh — ^- Then L{Sh) > {t — S)n. 
Proof. 

n n n 

L{Sh) > ^{Vh - Vi) = nyh-^yi> nyt - ^ \yi\ > nyu -6n> {t- 6)n. ^ 

1=1 1=1 i=l 

Lemma 3.6. Assume that \yh\ < i = 0.6. Letp £ S be a point in the middle strip Vm, with 
y-coordinate satisfying \y\ < 0.15. Then dmaxip) < 0.9605. 

Proof. It is straightforward to check that the maximum distance is attained for a point p 
on £2 with y-coordinate —0.15. The furthest point from p in the allowed region is c. Hence 



dmax(p) < \pc\ = \Jw'^ + (0.15 + t)2 = ^J{).& + 0.75^ < 0.9605. ■ 

We now distinguish the following four cases to complete our estimation of the approxi- 
mation ratio. 

Case I: Y^=i l^il — The algorithm output^ Sa or S\j. By Lemma [331 the approx- 
imation ratio is at least 



2L{Topt) " V 4 

Case 2: Y17=i IVil — ^''^ ™^ Vh ^ ^- The algorithm outputs S^. By Lemma [3.51 the 
approximation ratio is at least t — 5 = 0.55. 

Case 3: \yi\ < Sn and yh 1^ t and + n;, > (1 — z)n. The algorithm outputs 

Ea or El,. We only need the last inequality in estimating the length. By Lemma |3.4[ the 
approximation ratio is at least 

L{Ea) + L{Eb) ^ 1 + w 2n + na + n, {l + w){3-z) 1.6-2.52 

^ — > • > = = U.5U4. 

2L(Topt) ~ 4 2n - 8 8 

Case 4-' XliLi \yi\ ^ ^i^d yh 1^ t and Ua + Uf, < (1 — z)n. The algorithm outputs Sa 
or S}). There are at least zn = 0.48n points in the middle strip V^- Observe that at most 
n/3 points in Vm have \yi\ > 0.15; otherwise we would have 

n 

\yi\ > \y^\ > • ^ = 0.05n = 6n, 

i = l Vm 

a contradiction. It follows that at least 12n/25 — n/3 = lln/75 points in the middle strip 
have < 0.15. By Lemma [3.21 and Lemma [3.61 



64?T/ 11 77/ 

-^^(TbpT) < ^ + 0.9605 • — < 0.9943n. 

The approximation ratio is at least 

L(Sa) + L(Sb) n 

> > 0.502. 

2L(Topt) ~ 2 • 0.9943n " 

This completes the list of cases and thereby the proof of Theorem 11.21 

Remark. The example in Fig. [3] with n points (n even) equally spaced along a circle shows 
that the constant 0.502 measuring the approximation ratio achieved by our algorithm A3 



^Here and in other instances it is meant that the algorithm outputs a structure at least as long as these. 
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cannot be improved to anything larger than 2/7r. Indeed the lengths of the five structures 
computed by the algorithm are L{Sa) = L{Sb) = L{Sh) = L{Ea) = L{Ei,) = (1 — o{l))^n, 
while L{Topt) > L{Hopt) = (1 - o(l))n. 




Figure 3: The non-crossing structure Ea for an example with n = 16 points on the circle. The 
middle strip Vm is bounded by the two dashed vertical lines. 



4. The Hamiltonian cycle 

In this section we present the proof of Theorem ll.3t which is similar (including notation) 
to that of Theorem 11.11 The rotated coordinate system Ta, and the x-coordinates Xi with 
respect to this system are denoted in the same way. By relabeling the points assume that 
the optimal cycle is Qopt = Pi,P2,---,Pn (with the convention that Pn+i = Pi)- We 
approximate Qopt by constructing a non-crossing alternating path A on a subset of S, and 
then completing it to a non-crossing cycle using convex hull vertices. We need to observe 
that the alternating path A on the subset / of interior (non-hull) vertices of S produced 
by the algorithm of Abellanas et al. [I] is not good enough for this strategy: even though 
one endpoint of A (the first computed by the algorithm) is always on the convex hull of /, 
the other endpoint might be blocked by edges of A, so that A might not be extendible to a 
non-crossing Hamiltonian cycle (an example is shown in Fig. H}. Here, we give a stronger 
result that fits our purpose (for an even number of points). 

Lemma 4.1. Let S = RU B with with \R\ = \B\, be a linearly separable bipartition given 
by line £. Then S admits an alternating non-crossing spanning path A such that (1) the 
edges of A cross I at points ordered monotonically along i; and (2) the two endpoints of 
A are incident to the two distinct edges of the convex hull that connect R and B (the two 
red-blue bridges). Such a Hamiltonian path can be computed in 0(n log n) time. We refer 
to the underlying procedure as the two-endpoint path construction algorithm. 

Proof. We modify the algorithm of Abellanas et al. for path construction, so that the path 
is grown from the two endpoints and the two sub-paths merge "in the middle". Recall that 
S = RU B, and \R\ = \B\, thus \S\ is even. Let ribi and r262 be the top and bottom 
red-blue edges of the convex hull conv(S'), respectively, called top and bottom bridges; it 
is possible that ri = r2 or bi = 62 but not both. One endpoint of A is an endpoint of 
the top bridge, and the other endpoint of A is an endpoint of the bottom bridge, and they 
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Figure 4: A non-crossing alternating path obtained by the algorithm of Abellanas et al. For the 
purpose of cycle construction, the path is non-extendible from its 2nd endpoint, vertex 7. 

are chosen of opposite colors. Let A = {^1,62} 01 A = {61, arbitrarily, containing two 
endpoints of the path. At every step, recompute the top and bottom bridges of S" \ ^, 
and append either the red or the blue vertex of each bridge to A such that the appended 
edges cross the separating line £. In the last step, the convex hull of 5 \ A is a red-blue 
segment that merges the two sub-paths. The two new edges added simultaneously at each 
step cannot cross each other; and they cannot cross previous edges, since they are separated 
from them by the convex hull of S \A. Finally, they cannot extend the two sub-paths by 
the same point either, because |5| is even. ■ 

The next lemma follows from Lemma 2.1]; we will only need its corollary, Lemma [4.31 

Lemma 4.2. ([lO]). Let P = pi,P2, ■ ■ ■ ,Pn be a simple polygon (with the convention that 
Pn+i = Pi) and q be a point in the exterior of the convex hull of P, where P U {q} is in 
general position. Then q sees one edge PiPi+i of P. Such an edge can be found in 0{n) 
time. 

Lemma 4.3. Let P = pi,P2, ■ ■ ■ ,Pn be a simple polygon (with the convention that Pn+i = 
pi) and q be a point in the exterior of the convex hull of P, where P U {q} is in general 
position. Then the polygonal cycle P can be extended to include q so that P U {q} is still a 
simple polygon. More precisely, there exists i £ [n], so that Q = pi, . . . ,pi, q,pi^i, . . . ,Pn is 
a simple polygon. Moreover, L{Q) > L{P). The extension can be computed in 0{n) time. 

Proof. By Lemma 14.21 Q sees one edge PiPi+i of P. Replacing this edge of P by the two 
edges piq and qpi+i results in a simple polygon Q = pi, . . . ,pi, q,pi^i, . . . ,pn. By the triangle 
inequality, L{Q) > L{P). The extension can be computed in 0(n) time, as determined by 
the time needed to find a visible edge. ■ 

Note that the condition in the lemma that q lies in the exterior of the convex hull of 
P, is indeed necessary. Otherwise one cannot guarantee that q sees an edge of P. 

(i) Let S = S'uS" , where 5' is the set of convex hull vertices and S" is the set of interior 
points. Let S' = {pji,Pj2, ■ ■ ■ ^Pjh}- ^ = |5"|, m = 15"'!, thus n = h + m. Assume first 
for simplicity that m is even. An easy modification of the algorithm, explained below, is 
used for m odd. 
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Algorithm A4: 

Step 1. For all non-equivalent bisections of S" (i.e., for all balanced bipartitions of S"): 
1. Compute a non-crossing alternating path A by using the two-endpoint path construction 
algorithm (Lemma 14. ip . 2. Extend A to a cycle by connecting its endpoints to (one or two) 
convex hull vertices. 3. Further extend this cycle to include the remaining hull vertices, by 
repeated invocation of Lemma 14.31 

Step 2. Output the longest such cycle (containing all points of S). 

Observe that after Step 1.1, the two endpoints of the path are vertices of conv(S"'), 
hence they can be connected to hull vertices to make a cycle. If m is odd, then there is a 
point q £ S" on the line i. Use the two-endpoint path construction algorithm for S" \ {q}, 
and the same bisecting line i. If q is in the interior of conv(S"' \ {<?}), then extend the path 
with point q, using Lemma [2121 Otherwise, q sees the top or bottom bridge of conv(S'"\{g}), 
so the path can be extended by connecting q to the endpoint visible to q. The two endpoints 
of the extended path are on conv(5"), hence they can be connected to hull vertices to make 
a cycle, as in the case of even m. 
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